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of the renewal process, which have applications to bead-assay experiments performed 
both with processive molecular motors, such as myosin V and kinesin, and rotary 
motors, such as Fl-ATPase. 
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1. Introduction 

There has been a growing interest in the biology and biophysics communities in the 
study of the operation of motor proteins [Tj. Such biological agents are typically 
single enzymes that can act as thermodynamic engines, directly converting chemical 
energy into mechanical energy [2] through a chemical cycle which occurs at constant 
temperature 3J. Different classes of motor proteins are involved in biological processes 
such as cellular transport, cellular mitosis or muscle contraction, to cite a few examples 
of a plethora for which the presence of such agents is crucial. 

In this article, we concentrate on two particular classes of motor proteins, namely 
linear processive molecular motors and rotary molecular motors. Linear processive 
motor proteins, of which myosin V and kinesin jH Ej are important examples, catalyse 
the ATP hydrolysis reaction, ATP — > ADP+P, and use the energy so obtained to move 
along linear molecular tracks, carrying organelles or membrane patches with them and 
performing directed transport within cells with high efficiency [Oj. These organelles are 
too massive for molecular diffusion to move them efficiently in a crowded environment 
such as the cytoplasm within the time scales relevant for biological processes. The 
molecular tracks, which are composed of actin in the case of myosin V and of tubulin 
in the case of kinesin, have a polar character, i.e. such motors can move in only one 
direction. A different processive molecular motor, dynein, moves along the tubulin track 
in the direction opposite to kinesin. Yet another example of a linear processive molecular 
motor is the enzyme RNA-polymerase; see [Zj and references therein for a more complete 
discussion of the characteristics of this motor. Here it suffices to say that this enzyme 
moves along a DNA-strand, also using ATP hydrolysis as its energy source. However, 
the function of this enzyme is different from the motors described above, in that it does 
not perform molecular transport; instead it promotes the transcription of messenger- 
RNA from the underlying DNA-strand on which it moves. Another significant difference 
is that the motion of this motor occurs in an heterogeneous medium such as the DNA- 
strand, rather than in homogeneous media such as the actin or tubulin molecular tracks. 

The chemical cycle that results in the hydrolysis of one ATP molecule by the 
molecular motor is composed of several substeps, corresponding to changes in the 
internal conformation of the motor, accompanied by chemical reactions, e.g. ATP 
binding to the motor, ATP hydrolysis, ADP and P release, etc. Typically, one ATP 
hydrolysis is necessary for the motor to advance one step along the track, docking to 
the next available site in it, the size of such steps being around 35 nm in the case of 
myosin V and 8 nm in the case of kinesin [HI El QUI E] • It should also be noted that 
the use of the word 'processive' to classify a molecular motor is limited to motors whose 
dwell time on the molecular track is much larger than the time for a complete chemical 
cycle, i.e. the typical time for a single step along the track. 

Another class of molecular motors to which one may also apply some of the results 
presented in this paper are the so-called rotary molecular motors, of which Fl-ATPase 
is perhaps the most important example j!2j . This motor can also catalyse the ATP 
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hydrolysis reaction and use the energy so obtained to generate the rotary motion of a 
shaft-like mechanism. The mechano-chemical cycle of such a motor is composed of three 
rotations of 120 degrees, determined by the symmetry of the molecule in question, each 
rotation being coupled to the hydrolysis of one ATP molecule, as described by the so- 
called "binding change mechanism" ^3] ■ If such a motor is coupled to the proton-driven 
Fo motor, the assembly of these two motors can either work as a proton pump or use a 
proton gradient to synthetise ATP, hence its name, ATPsynthase. 

A great wealth of knowledge about the internal chemical cycle of linear processive 
molecular motors has been acquired through bead-assay experiments [HJ [T3J |UJ HEl HI 
El El EE]- In such an experiment, performed in a fluid medium, a molecular motor 
is coupled to a dielectric bead of micrometric size, the size of such a sphere being 
nevertheless much larger than that of the motor protein. The bead can be manipulated 
through the use of optical or magnetic tweezers. In certain situations, one can, using 
an optical trap, perform experiments in which the molecular motor, moving along its 
track, is subjected to a known constant force of several pN. The study of the motion for 
different values of the applied force allows one to determine the force-velocity relation 
characteristic of a given motor. At a given value of the force, known as the stalling 
force, the motor usually decouples from the track (such decoupling can even occur in 
the absence of force, but it is a much rarer event in that case ^H]) or it may move 
backwards along the track. Note that such experiments may be performed at distinct 
values of the applied load and also at different ATP concentrations. 

More recently, Cappello et al. l2"T] have considered an experimental apparatus 
containing a bead-motor assay that moves through the interference fringes of an 
evanescent light-wave that exists in the proximity of a microscope glass plate. The 
force which the electric field of such wave exerts on the bead is negligible, i.e. the 
experiment is performed in a zero-load condition. However, the bead still scatters 
photons of the evanescent field and the observation of such events may be used to 
track the bead's position with high spatial and time resolutions. In particular, such a 
method allows for measurements with time resolutions of [is (MHz). In normal bead- 
assays, the feedback mechanism used to keep the bead under constant force typically 
limits the time resolution to a few milliseconds (kHz). Note that one may perform these 
measurements at different ATP concentrations, but one is always limited to work at 
zero external load, which constitutes both a strength (because it permits higher time 
resolutions) and a weakness of the method (because it limits the parameter range in 
which the system can be studied). 

In essence, a clear qualitative picture has emerged from the different types of bead- 
assay experiments, namely that a few of the substeps which compose the chemical 
cycle of a molecular motor (known in the literature as the rate-limiting steps) are of 
particularly long duration compared to the remaining substeps, and that the statistics 
of the motion of the molecular motor are, for the time resolutions available, essentially 
determined by the duration of such rate-limiting steps. Thus, the experimental study 
of the motion of linear processive molecular motors can provide valuable information 
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concerning the nature and the duration of the rate-limiting steps under different load 
conditions or under different concentrations of ATP ^H] ■ There is experimental evidence 
that two rate-limiting steps are sufficient to describe the chemical cycle of myosin V, 
namely ATP binding (by lowering enough the concentration of ATP it is always possible 
to make ATP binding to the motor a rate-limiting step) and ADP release jlj, whereas 
three or four rate-limiting steps may be needed to describe the chemical cycle of kinesin, 
depending on the ATP concentration. It is not yet well understood to which chemical 
substeps these rates correspond [T%] . 

The rotary movement of the 7 sub-unit of the Fl-ATPase motor can be visualised 
by coupling such a motor to a microprobe, which can be either a fluorescent actin 
filament, a gold or polysterene bead, or a single fluorescent dye [T21 l22j. The rotation of 
the probe can be recorded with a microscope and a CCD camera. It follows from such 
studies that the 360-degree rotation can be decomposed into three 120-degree substeps, 
each of which is in turn composed of two rate-limiting steps with approximately the 
same duration. These steps involve rotations of 90 and 30 degrees, respectively [23J. 
One should note that of the three different types of microprobe experiments referred, 
the use of a single fluorescent dye seems to be the most promising method of visualisation 
J21I22I, as it does not involve a perturbation of the motion of the molecule. In the other 
cases, the rotation of Fl-ATPase is hindered by the large frictional coupling between 
the microprobe (actin filament or gold or polysterene bead) and the surrounding fluid. 

The theoretical tool that we will be applying to the analysis of such experiments is 
the concept of renewal process. Such processes are ubiquitous in physics. These models 
have been successfully applied to describe the motion of processive molecular motors 
such as myosin V and kinesin [HI [201 HI] , the statistics of detection of quantum particles 
[2U 123126] and persistence phenomena in kinetic Ising models [2Z1I2E]- Also, the results 
that relate to persistent phenomena have applications to studies of the volatility of 
financial markets [20] • The above list is not exhaustive. 

Loosely speaking, a renewal process is a counting process where unit increments 
occur at random times. We will consider only independently distributed renewal 
processes, i.e. processes in which the probability distribution for the occurrence of 
the next increment (also known as the waiting-time distribution of the renewal process) 
depends only on the time elapsed since the occurrence of the last increment and not on 
the previous history of the counting process. 

If one wishes to be more specific, one can say that the purpose of this article 
is to enumerate and classify a series of bead-assay experiments performed with linear 
processive molecular motors or with rotary molecular motors. These experiments are 
characterised by the fact that one can apply the simple (and exact) results obtained 
from the calculation of multiple-time correlation functions of renewal processes with an 
arbitrary waiting-time distribution to the analysis of the data obtained. In order to 
apply such simple models to linear processive motors, one represents a forward step of a 
molecular motor by an increment of the renewal process, whose waiting-time distribution 
is dependent on the number of rate-limiting steps of the motor. Such an approximation 
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is justified if the molecular motor performs backward steps only infrequently within 
the time-window of observation, i.e. provided one is not working too close to the stall 
force or at ATP concentrations which are too low, when the rates for a forward or a 
backward step to occur become comparable and one needs to use a large time-window of 
observation. The application of these results to the dynamics of rotary motors involves 
the mapping of a rotation of the motor by a given angle to an increment of the renewal 
process, but such a mapping is more subtle (see below for details). Such a mapping 
can be justified provided that the motor performs backward rotations only infrequently 
within the time-window of observation, which is the case if one is working at high ATP 
concentrations. 

The measurement of single-time correlation functions of the number of increments 
N(T), which have occurred in a renewal process until a given time T, is a common 
practice in the context of the experimental study of linear processive molecular motors. 
The measurement of these quantities was first undertaken by Svoboda and coworkers 
in their experiments performed with the linear processive molecular motor kinesin |14j . 
These authors have considered the behaviour of the first and second moments of N(T), 
i.e. (N(T)) and (N 2 (T)), where the averaging is taken over different realisations of 
the experiment. Such correlation functions contain information concerning both the 
number of rate-limiting steps in a chemical cycle and the characteristic rates pertaining 
to such steps. However, such single-time correlation functions do not fully characterise 
the motor's chemical cycle, in particular in the case of motors whose cycles are composed 
of many rate-limiting steps, like kinesin. The measurement of multiple-time correlation 
functions can provide additional valuable information in such a case fK)\ I2T] . 

We employ a method based on the use of the probability-generating functional 
PHI EH] to compute multiple-time correlation functions of a renewal process. For 
simplicity, and due to its experimental relevance, we explicitly compute the mean- 
square deviation of the number of increments that occur between time t 2 and a later 
time t\. This correlation function is mathematically defined as {(N(t\) — N(t 2 )) 2 ) — 
(N(ti) — N(t2)) 2 , where the times t 2 and £1 are large compared to the typical time of 
a single chemical cycle. We show that this function contains additional information 
concerning the rate-limiting steps of the chemical process, information which cannot 
be extracted from single-time correlation functions. We also briefly indicate how the 
computation of higher-order correlation functions can be performed. Furthermore, we 
explicitly compute the density-density correlation function for a single motor, which was 
considered in the experiments of Cappello and co-workers Interestingly enough, 

one can show that in a certain limiting case, such a function is identical to the spin-spin 
correlation function considered by Godreche and Luck |2Zj in their study of persistence 
phenomena in kinetic Ising models. Furthermore, the class of spin models obtained in 
this particular limit (to which the model studied by Godreche and Luck belongs) may 
be relevant for the experimental study of rotary molecular motors, and hence it is also 
discussed here. 

The structure of this paper is as follows: in the next section, we present the general 
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results obtained and relate such results to the relevant experiments. In section we 
provide a mathematical introduction to renewal processes, following [3U| and introduce 
the probability-generating functional of the renewal process, together with some related 
quantities whose usefullness will become apparent in section HJ In section HJ we provide 
a general derivation of the results presented in section concerning multiple-time 
correlation functions of a renewal process, including the mean-square deviation and 
the density-density correlation functions mentioned above, whose expressions we will 
explicitly compute in the asymptotic regime of large times. A reader whose primary 
interest is not mathematics may skip sections H3 and 0] without loss of continuity with 
respect to the remainder of the paper. Finally, in sectional we present our conclusions 
and provide a brief outlook of the experimental work that we believe can be carried 
in this field using the results that we have derived. In the appendices, we discuss two 
simple examples of waiting-time distributions for which the quantities discussed in the 
main text can be computed outside the asymptotic regime. 

2. General discussion 

The most elementary quantity one can measure in a bead-assay experiment is the average 
displacement of a motor (x(T)) up to time T over many runs. The simplest model for 
such a motor is that x(T) = N(T) d, where N(T) is an integer variable, updated at 
random intervals with a given distribution f(() (a renewal process), and d is the motor's 
step size. In the limit of large time, Blackwell's renewal theorem [32] guaranties that 
(x(T)) = T dj (t) asymptotically, where (r) = d((f(Q is the average waiting time 
of the distribution and corresponds to the mean duration time (or turnover time) of a 
single chemical cycle. 

The simplest choice one can take for f(Q is the exponential distribution f(Q = 
e~ c -l T jr. In this case, (r) = r. In their discussion of bead-motor assay experiments 
performed on the processive motor kinesin, Svoboda, Mitra and Block ^3] have 
examined the more complicated situation where the waiting-time distribution f(() is 
given by the convolution of a finite number A4 of simple Poisson processes. Each of these 
processes, which occur in series with typical times r 1; • ■ ■ , r M , is supposed to represent 
a rate-limiting step of the motor's chemical cycle. In such a case, the average waiting 
time is given by 

M 

<r> = (1) 

i=l 

Such a model is in good agreement with experiments where one measures the distance 
travelled by the molecular motor alone. If one wishes to resolve the chemical 
cycle substeps and its associated pathways, e.g. through the use of cryoelectron 
microscopy p£H or by using bead-motor assays where lateral or forward loads 
are applied to the motor ^Hj, one needs to make use of more involved models 
[SH IS3 El ISEl EZl EH! EHJ IHH EU in order to interpret such experiments. 
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Svoboda et al. [U] have also introduced the concept of the randomness coefficient 
r, characteristic of a single molecular motor. Given that the quantity (x 2 (T)) — (x(T)) 2 
is the mean square deviation of the distance travelled by the motor, r is defined as 



The randomness coefficient, being the ratio of the mean-square deviation of distance 
travelled by the motor to the average distance travelled itself, represents a measure of 
the deviation of the motors stepping from a deterministic motion, which would occur 
if r = 0. Also, note that r is chosen to be dimensionless |42| . For the simple Poisson 
process, (x 2 (T)) — (x(T)) 2 = Td 2 /r, i.e. r — 1. For a renewal process composed of M. 
rate-limiting steps [13], one can show ^1] that the randomness parameter is given by 



It follows from (J5j) that if A4 > 1, r < 1. For a renewal process composed of Ai rate- 
limiting steps, r > the equality being obtained when all rates are equal. In the 
limit of an infinite number of substeps whose characteristic time tends to zero, r = 
and the motor performs a deterministic motion. 

The measurement of the randomness parameter, which is robust against thermal 
noise ^1] or the influence of the initial conditions (i.e. of the experimental set up), is a 
powerful experimental tool that can be used to rule out a proposed chemical cycle, if such 
a cycle contains too small a number of rate-limiting steps. For example, a measurement 
of r < 1/2 indicates that at least three rate-limiting steps are needed to describe the 
motor's chemical cycle [4*4*] . 

Since the chemical cycle of myosin V appears to be composed of only two rate- 
limiting steps, the joint measurement of r and of the average distance travelled by the 
motor at large times is sufficient to determine the value of the typical times T\ and 
T2 associated to each substep. However, if one is considering the experimental study 
of motors with more than two rate-limiting steps in their cycle, such as kinesin, the 
measurement of these two quantities is not equivalent to the complete determination of 
the characteristic times . . . , r M (M. > 2). 

The statistical analysis of step duration using such bead-motors assays, which allows 
for the direct extraction of the waiting-time distribution /(C), was also considered by 
several authors [HE]. By fitting f(Q with an appropriate convolution of exponential 
functions, one can obtain the values of all the characteristic times t±, . . . , tm for arbitrary 
A4. However, such a technique requires the substitution of a signal with a rich structure 
(the individual trajectories of the motor as a function of time) by a function in steps, 
which leads to an effective truncation of the data, because possible substeps within 
the chemical cycle are erased. Furthermore, such a technique requires a resolution of 
individual steps, which is not always feasible. 

One is thus led to consider the information provided by multiple-time correlation 
functions. The simplest quantity one can consider is the average number of steps given 




(2) 



r = 




(3) 
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by the motor between times t 2 and ti, i.e. the correlation function (N(ti) — N(t 2 )), 
which for large times t±, t 2 — > oo behaves as (ti — t 2 ) / (t), and which obviously does not 
carry any new information. Its measurement yields, as above, the value of the turnover 
time (r) of the chemical cycle. 

The same cannot be said of the mean-square deviation of this quantity, i.e. the 
connected correlation function ((Nfa) - N(t 2 )) 2 ) conn = ((N(t x ) - N(t 2 )) 2 ) - (Nfa) - 
N(t 2 )) 2 . If one were to take t 2 = 0, this correlation function would reduce to the mean 
square deviation of the distance travelled by the motor up to time t\, introduced above, 
and from which one can extract the randomness parameter, as defined in (J2J). However, 
there are good reasons to consider instead the opposite limit in which both ti,t 2 — ► oo, 
with t = t\ — t 2 kept finite, as in this limit the correlation function becomes dependent 
only on t, i.e. one recovers a form of time-translation invariance (we will write the 
correlation function in this limit as C(t) = \im tli t 2 ^oo((N (t\) — N (t 2 )) 2 ) conn ) . Such a limit 
is of experimental relevance, since such a correlation function becomes independent of 
the initial conditions, which can change from one experimental realisation to the other. 
Furthermore, its measurement contains additional information concerning the waiting- 
time distribution, rather than just the value of the randomness parameter, as we will 
now discuss. Finally, such a measurement, just like the measurement of the average 
number of steps (N(ti) — N(t 2 )), requires only the ensemble averaging of the data and 
is not plagued by the limitations one encounters if one tries to measure f(() directly, as 
described above. 

In our treatment, we will consider a renewal process with an arbitrary waiting-time 
distribution, rather than the special choice made by Svoboda et al., since the study of 
the general case does not involve a more complicated analysis. However, we will always 
indicate the results for this particular waiting-time distribution, given its experimental 
importance. 

It will be shown in section 0] that this correlation function displays the following 
behaviour in the limit of small and large time differences t = t\ — t 2 , 

\ j\ if t < (r) 

{-fa t + c lf t »( r ) 

with r = ( r2 ]~fi and C = - ^X, where (r 2 ) = / °° rfCCV(C) and (r 3 ) = 

(r) 2(r) s(r) 

Jo°° d(( 3 f((). Please note that the result just quoted for r, the randomness parameter, 
is a generalisation of equation (jSJ) to the case of general waiting-time distribution f(() 
[27] . The constant C can be written in terms of r and of the connected third- moment 

of the waiting-time distribution, 5 



((r-(r)) 3 ) 



1 1 2 

C = - + -r 2 --6. (5) 
6 2 3 v ; 

In the case in which f(() is given by the convolution of A4 Poisson processes, as 

considered by Svoboda et al., one can show (see section HJ) that r reduces to equation 
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(J3J). In this case, 5 



V^ 1 t3 



-. The constant C is exactly zero for a simple Poisson 



process. 

As we will show in |Appendix A[ one can compute C{t) for arbitrary values of t, 
rather than just the asymptotic limits given by (0J), for processes whose waiting-time 
distribution is given by the convolution of two or three Poisson processes, i.e. if = 2,3 
(the case M. = 1 is trivial, see section EJ). Such computations are explicitly performed 
because of the importance of such processes for the experimental study of myosin V and 
kinesin. We plot the results of these calculations below, in figures [T] and El respectively, 
with the choices r± — t% — 1/2 for the case of a renewal process composed of two rate- 
limiting steps and T\ = T2 = r% = 1/3 for the case of a renewal process composed of 
three rate-limiting steps. In both cases, (r) = 1 (in arbitrary units), with r = 1/2 in 
the first case and r = 1/3 in the second case. It is seen that the two functions have the 
correct asymptotic limits at short and large times, as given by equation (@J). 



Correlation 
2 




1.5 2 



Figure 1. Mean-square deviation of the number of steps in the regime of large times 
(red curve), for a motor with two rate- limiting steps. Also shown are the linear regimes 
at small and large time t. 




Figure 2. Mean-square deviation of the number of steps in the regime of large times 
(red curve), for a motor with three rate- limiting steps. Also shown are the linear 
regimes at small and large time t. 



As stated above, the measurement of the average distance travelled by the motor 
permits one to determine the turnover time (r). Once this quantity is known, one 
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can use the measurement of C(t) at large times to determine both the randomness 
coefficient r (which is given by the slope of the straight line multipled by (r)), as well 
as the constant C which is related to the second and third moment of /(£) j3H] . For a 
motor whose chemical cycle is composed of two rate-limiting steps, the constant C is 
given by C = p^^i ■ The measurement of C does not provide any new information 
in this case, but it may provide a way to check or to improve the results obtained 
from the joint measurement of the average distance travelled by the motor and of the 
randomness parameter r. On the other hand, in the case of a motor whose chemical cycle 
is composed of three rate-limiting steps, the measurement of the total distance travelled 
by the motor, which permits one to determine the total turnover time (r) = T\ + t 2 + T3, 
of r and of C — ( Ti+T ^ +r3 ^ (rf r| + rfr 2 r 3 + Tirf r 3 + rf t| + t\t 2 t1 + rf rf ) can be used 
to determine the three time-constants n, t 2 and T3, through a fit of the experimental 
results. Such a measurement may be particularly useful in experiments with kinesin and 
for concentrations of ATP for which ATP binding is not a rate-limiting step, because in 
such a case the chemical cycle of kinesin appears to be composed of three rate-limiting 
steps [THj . 

Another example of a multiple-time correlation function which can be measured 
experimentally is the density-density correlator of an ensemble of independent molecular 
motors, defined as S(q,t 1 ,t 2 ) = ( e~ tqd ^ N ^ tl ^~ N ^ 2 ^ ), where q has the dimensions of a 
wave-vector and d is the motor's step size. This quantity is the Fourier transform of the 
probability for the motor to move by a distance n x d between time t 2 and t±. It was 
directly measured in the experiments of Cappello and coworkers jlHj, where the bead- 
motor assay moves through the interference fringes of an evanescent wave. In such an 
experiment, q = 2tt/A, where A is the period of the mask used to create the interference 
pattern [20J |2T]. In the long time limit t 2 ,t 1 — > 00, S(q,ti,t 2 ) becomes solely dependent 
on the time difference t = t\ — 1 2 , in which case one simply writes S(q, t). One can show 
(see sectionHJ) that the Laplace transform of such a function, S(q, s) = / °° dt e~ st S(q, t), 
is given, for a general distribution f(() with a finite average time (r), by 

^.l^fH^-""-^), (6) 

where f(s) is the Laplace transform of /(C)- Note that in the experiments of Cappello 
and co-workers, this function was taken to be a simple exponential distribution. The 
equation above generalises their result to the case of a general waiting-time distribution. 

One can relate the Fourier transform of S(q,t), S(q,u) = dt e luJt S(q,t) to its 
Laplace transform by S(q,u) = 2ReS(q,s = —iuj). The Fourier transform is more 
amenable to computation from the measured data than the Laplace transform, or can 
otherwise be directly measured. One can show from such a relation and from equation 
(P) that S(q, uj) is given by 

V ; u 2 (T)\l-e-^ d f{-iu)\ 2 V ; 
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If /(£) is a simple exponential distribution, one obtains from (J7J, 

{q ' } (cut - sin(grf)) 2 + (1 - cos(gd)) 2 ' W 
which becomes the result for asymetric diffusion in the continuum limit qd <C 1, i.e. one 
obtains in this limit 

where we have introduced the motor's velocity v = d/r and the motor's diffusion 
constant D = d 2 /2r, valid for a Poisson process, since equation (jHJ) was derived from 
(JSJ). For a renewal process with waiting-time distribution /(C), one can still approximate 
S(q,u) by © if the Laplace transform f(s) is analytic at s = 0, i.e. if all moments 
of the distribution /(£) exist. In such a case, v = d/(r) and D = d 2 r/2(r). Such 
approximation is valid at low-frequencies compared to the total turnover rate and at 
low-momenta compared to the inverse step size, i.e. if u{t) <C 1 and qd <C 1. In 
other words, the motor behaves in such a limit as a Brownian particle, characterized 
by the two parameters v and D 21], as one would expect. In |Appendix B[ we will 
explicitly compute S(q,u) outside this asymptotic region for the cases in which /(£) 
is given by the convolution of two or three Poisson processes, which corresponds to a 
chemical cycle composed of two or three rate-limiting steps. Again, such calculations 
are carried through because of the experimental significance of these two cases for the 
study of myosin V and of kinesin, since a measurement of S(q, t) and the subsequent 
computation of S(q, u) would permit the extraction of the relevant time constants by 
fitting the measured S(q,u) with the corresponding expression, appropriate for the 
given number of rate-limiting steps, as given in |Appendix B| 47]. Note that in order to 
perform such a fitting, one should first have a qualitative understanding of the nature 
of the chemical cycle, in particular of the number of rate-limiting steps. 

Interestingly enough, one may also use the results obtained above for S(q, s) to 
interpret experiments performed with rotary molecular motors, such as Fl-ATPase 
[T2*l IT3] . In this case, one assimilates the rotary motion of the probe attached to the 
motor to the motion of classical spin in two dimensions (an XY model), with unit length 
and Q internal states, where Q is an integer, such that the spin rotates around the z-axis 
in one direction by an angle equal to 2tt/Q. Such a model may be appropriately referred 
to as a 'random-clock'. 

In more mathematical terms, one defines the 2-dimensional vector variable oV = 
(cos(27rA r (T)/Q), sin(27rA r (T)/Q)), where N(T) is given by a renewal process with an 
arbitrary waiting-time distribution /(C)- Such a vector rotates in the 2-dimensional 
circle by an angle equal to 2tt/Q, with the waiting-time distribution of the renewal 
process /(C)- After Q rotations, the spin returns to its original state. It can 
be easily seen that the spin-spin correlation function of such a model is given by 
{&t\ " &t 2 ) = ReS(q = 2ii/Q,ti,t2), where S(q,ti,t 2 ) is the density-density correlator 
defined above for the renewal process with waiting-time distribution f(() and d = 1, 
since q is here an angular variable. In the long-time limit, ti,t2 - * oo, where this 
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function is solely dependent on t = t\ — t^, one can read the Laplace transform of such 
a quantity from equation (JSj). This transform is given by 




(10) 



(r)s(l-e- 27rl /Q/( s )) / 

If we were to take Q — 2, the model would become that of an Ising spin variable and 
equation (fTUj) reduces to the result obtained by Godreche and Luck [27]. In effect, the 
computation of S(q, s), as given by (jHJ), can also be obtained from their results. We will 
use a different method to derive such an expression. 

In the case of Fl-ATPase, the symmetry of the motor molecule implies that Q = 3, 
i.e. the motor rotates by an angle of 120 degrees and requires three such rotations to 
complete its cycle. If one were to assume that the waiting-time distribution for each of 
these 120 degrees rotations was given by a simple Poisson process, it becomes trivial to 
invert the Laplace transform given in equation (fTU|) and one obtains (with t\ > t 2 ) 

i.e., the spin-spin correlation function, despite the irreversible character of the underlying 
spin model, displays an oscillatory, albeit overdamped, behaviour. In fact, such 
behaviour is always present if Q > 2. Its origin is trivial, being traceable to the 
circular character of the motion. Nevertheless, the study of the motion may provide 
valuable information concerning the chemical cycle. In equation (jHJ, such information 
is encoded in the turnover time r. 

In reality, the 120-degree rotation performed by Fl-ATPase is composed of two 
rate-limiting substeps, the first substep corresponding to a 90-degree rotation and the 
second substep to a 30-degree rotation This implies that f(Q is given by the 

convolution of two Poisson processes. In such a case, it is still possible to invert the 
Laplace transform given by equation ([lip, but the resulting expression for the spin- 
spin correlation function in the time domain is rather cumbersome and it is preferable 
to work with its Fourier transform instead. However, there is one particular instance 
in which the spin-spin correlation function in the time domain acquires a particularly 
simple form, namely when the time-constants t% and t 2 are equal, i.e. T\ = t 2 = r/2. 
One obtains for the spin-spin correlation function in this case, the result (with t% > t 2 ) 

,-, - v fV3(tl-* 2 )\ /3 _(ti-ta) 1 _3( tl -* 2 ) \ 

(o- tl ■ Ofc) = cos I J l^-e r +_ e r j, (12) 

whose form is clearly different from the one given in equation (jlljl . as the period of 
the oscillations is now shorter than its decay-time, i.e. the oscillations are no longer 
overdamped. This property may actually permit the experimental observation of such 
oscillations. We leave the derivation of equation (fT2*|) to |Appendix C We plot below, 
in figure H3 the two functions given by (fTTj) and ()12j) for comparison, in terms of the 
variable t = t\ — t 2 . 

The derivation of the expression of the Fourier transform of the spin-spin correlation 
function, for a general choice of the two time constants T\ and t 2 , is also left to 
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Spin-spin correlation 




-0.2 



Figure 3. Spin-spin correlation functions as given by equations l|llfl (blue curve) and 
(|12H (red curve). It is seen that the second function decays more slowly than the first 
and it also oscillates more rapidly. 



Appendix C Below, in figure 0] we plot such a function for the three distinct choices 
Ti = 1, r 2 = 0, which is just the Fourier transform of equation (fTTj). T\ = 4/5, r 2 = 1/5, 
and 7~i = r 2 = 1/2, which is the Fourier transform of equation (fT2*|) . The choice of units 
is such that the total turnover time r = T\ + r 2 = 1. Note that the difference between 
the two extreme cases T\ — 1, T2 = and t% = r 2 = 1/2 is more pronounced if the 
functions are plotted in Fourier space (figure EJ) rather than in time-domain (figure E]). 
The experimental situation is closer to the third case considered, i.e. the two time- 
constants are approximately equal to 1ms [12]. Here again, the measurement of the 
spin-spin correlation, followed by a Fourier transformation of the signal, is the preferred 
method for the extraction of the chemical cycle's time constants, since one can fit the 
experimentally measured signal with the result of equation (jC.3|) . which then yields T\ 
and r 2 . 

As we have mentioned above, the measurement of the spin-spin correlation function 
should ideally be performed using a fluorescent dye, since such a probe does not affect 
the rotation of Fl-ATPase. The visualisation of the single fluorophore, coupled to the 
rotating 7 sub-unit of Fl-ATPase, is performed by excitation of the fluorophore with 
circular polarized light and by measuring the intensities, H and V, of the light emitted by 
the molecule along two perpendicular directions of polarisation in a dual-view apparatus 
[22] • The geometry of the experimental setup is such that the polarisation of the observed 
light is parallel to the dipolar moment of the fluorophore, and is therefore determined by 
the instantaneous angle of rotation 6(T) = ^ N(T), of the 7 sub-unit with respect to the 
fixed a(3 sub-units [22J, where we take N{T) to be given by a renewal process, as above. 
Hence, one concludes that H(T) oc cos 2 (^V(T) + 0), V(T) oc sin 2 (^V(T) + 0), where 
(f) is the angle between the horizontal polarisation direction and the nearest polarisation 
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Figure 4. Fourier transform of the spin-spin correlation function for a renewal process 
composed of two Poisson processes. The three cases shown correspond to the choices, 
ti = 1, T2 = (blue curve), t\ = 4/5, T2 = 1/5 (green curve) and t\ = t-i = 1/2 (red 
curve). Note that the area under the curve is the same for all plots. 



direction along which light is emitted, with <ft = 18° in the experiments of Adachi and 
co-workers [221 - 

With such definitions in hand, we define h(t 2 ) = H(tf)+v(t 2 ) ' = 

- L \J H(t2)+v(t 2 ) ' sign in these definitions being determined by the value of the 
instantaneous polarisation P(t 2 ) = h 2 (t 2 ) — v 2 (t 2 ) [48] . Using the same definition for 
h(ti), v(ti) at a later time t 1 , one concludes that 

(a tl ■ a t2 ) = (cos^Nih) - N(t 2 ))/3}) = {h{h)h{t 2 )) + (v{ti)v{h)), (13) 

an equality that shows that the spin-spin correlation function can be related to the 
average value of a directly measurable quantity. One should note that the experiments 
of Adachi et al. [22] were performed at low ATP concentrations where ATP binding is 
a rate-limiting step, with typical binding times of the order of seconds. If one wishes 
to measure the characteristic times T\ and r 2 , one needs to perform the experiment at 
much higher ATP concentrations and with a time-resolution which is three orders of 
magnitude higher than in the experiment of Adachi and co-workers (n ~ t 2 ~ 1 ms). 
Achieving such a resolution may constitute an insurmountable practical problem, even 
if the method theoretically permits an unlimited time-resolution [4~§] . 



3. Definition of the renewal process and associated quantities 

A renewal process is defined by a random variable N(u) that is incremented by one unit 
at random times u = T\ < T 2 < . . . < where, in the simplest case and to which 
we will attain, the length of the intervals Ti,T 2 — Ti, . . . , TV — TV-i is taken from a 



Renewal processes and fluctuation analysis of molecular motor stepping 



15 



given waiting time distribution /(£) such that /(£) is positive and J °° d( f(() = 1 (one 
assumes that the process starts at u = with N(0) = 0). Mathematically, this can be 
written as 

N(u) = n if T n < u < T n+1 , (14) 

with Tq = 0. If one defines F k (T) to be the probability that at least k increments have 
occurred up to time T irrespective of whether more increments have occurred or not, 
then it is easy to show [32J that F k (T) is given by 

F (T) = 1 (15) 

F k (T) = [ T d(F k ^(T-()f(() k>0. (16) 
Jo 

One may equally regard the renewal process as a totally asymetric random-walk of a 
structureless particle, in which case one speaks of the probability of observing at least 
k forward steps of such a particle up to time T. As discussed above, this particle is 
supposed to represent a molecular motor, with the waiting time-distribution /(£) being 
determined by the internal chemical cycle of the motor. 

In particular, it follows from this definition that Fi(T) = Jq d( f(Q is the 
cumulative probability to observe an increment of N(u) in the time interval [0, T[. The 
probability that exactly k increments have occurred up to time T is then given by 

P k (T) = F k (T) - F k+1 (T) . (17) 

The average number of increments of N(u) (or the average number of steps of the 
particle) up to time T, which we designate by m(T), is given by m(T) = (N(T)) = 
Y^kLo kP k (T). Using equation (fTTjl for P k (T) and expressing F k+ i(T) in terms of F k {u) 
at earlier times through equation (|16jh it follows that m(T) necessarily obeys the integral 
equation 

m(T) = F X (T) + [ T d(m(T - C) /(C) , (18) 
Jo 

which is known in the literature as the 'renewal equation', m(T) being called the 'renewal 
function'. For reasons that will become clear below, we prefer to work with the derivative 
of m(T) and we write l(T) = m'(T), which we will also call the renewal function. 
Differentiating the above equation with respect to T and given that m(0) = 0, one 
obtains 

J(T) = /(T)+ [ T d(l(T-()f((), (19) 
Jo 

i.e. the equation has the same form as (|18|). but the non-homogeneous term is now 
given by f(T) rather than by F\(T). From the above equation, it also follows that 
1(0) = /(0). If one considers a Poisson process, where f(() = e~^ T /r, it is trivial 
to verify that l(T) = 1/r is the solution of equation (JTHJ) and one has m(T) = T/r, 
i.e. the average number of steps of the walker increases linearly with time. For more 
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general forms of the waiting-time distribution, the renewal equation can still be solved 
by Laplace transformation, and one obtains 

If s ) = f ( S ) , (20) 
1 - f(s) 

where l(s) and f(s) are, respectively, the Laplace transforms of l(T) and of f(T). In 
order to obtain l(T) from this expression, one needs to invert the Laplace transform, an 
operation one can only perform in a limited number of cases. Nevertheless, Blackwell's 
renewal theorem [32] asserts that provided that (r) = J °° d((f(() is finite, then 
l(T) = asymptotically at large time T. 

Having defined the random process that will be the object of our study, we now 
proceed to define the quantities we wish to compute. We will first consider the restriction 
of the renewal process given above such that exactly Af increments have occurred 
between time u = and time u = T, i.e. the times of the increments are such that 
< T% < T 2 < . . . < TV < T < TV+i and we will consider correlation functions of the 
random variable |3~T] 

Af 

J{u)=Y,9{u-Ti), (21) 
i=i 

where < u < T and where g(u) is a function with at most a finite number of 
discontinuities. If Af = 0, we define J{u) = for all < u < T. Note that 
N{u) = E^i 0{u - Ti), where 6{u) = if u < 0, 9(u) = 1 if u > 0, is of this form, and 
therefore the number of increments N(u) (or the number of steps which the particle has 
given) corresponds to a particular case of ()21|) in which one takes g{u) = 8(u). Thus, 
we wish to determine the value of the correlation functions (J{ti) J{t<i) . . . ff(t m ))j^, 
where m > 1 and where the averaging is over the distribution of the times T\ to T^-. 
For simplicity, we will restrict ourselves to m = 1 and m = 2 in the particular examples 
discussed, though the formalism is valid for arbitrary m. 

If one defines the probability-generating functional Zj^\j] as 

ZM = ^xp £ dxj(x) J(x)^j J , (22) 

where the averaging is performed over the increment times T%, . . . , TV, then it follows 
that the computation of Z^[j] is equivalent to the computation of all correlation 
functions, which can be obtained by functionally differentiating Zj^[j\ with respect to 
j. Note that we choose the normalisation of Z^\j\ to be given by 

Z M [j = 0] = P N (T) , (23) 

where P^/(T) is the probability of occurrence of exactly Af increments up to time T, 
introduced above. 

One can now generalise the concept of probability-generating functional to 
situations in which the total number of increments is not fixed by introducing the 
(grand-canonical) functional [3*T] 

oo 

m = e z m ■ ( 24 ) 

A/=0 
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The functional Z[j] can again be functionally differentiated with respect to j in order 
to obtain the relevant correlation function. Note that if one takes j = 0, one obtains 

oo 

Z[j = 0]= J2 Pn(T) = 1, (25) 

Af=0 

i.e. the probability-generating functional is equal to one due to normalisation of the 
total probability. This result justifies the normalisation chosen above for Zj^\j\. In the 
next section, we will develop a technique which will allow us to compute the correlators 
of J{u) order by order. 

4. Integral equation for the probability- generating functional and 
associated multiple-time correlation functions 

We first consider the case in which exactly Af > increments have occurred up to time 
T. The exclusive probability density (also known as Janoussi local density) EH] for 
the increments to occur at times T\ < T 2 < . . . < < T is given by 

pxiTu T 2 , . . . , T N ) = P (T - T M )f{T M - T^-i) ■ ■ ■ f(T 2 - T\)/ (T1) , (26) 

where Pq(T — TV) = 1 — J T ~ Tm d( f(() is the probability that no increment has ocurred 
between TV an d T, as given by equation (|T7jl (with k = 0). The probability-generating 
functional Zj^[j] is given in terms of this quantity by 

ZM= f ^ r "dTjV-i..- r dT lPM (T l ,T 2 ,...,T J ^)e i Io 
Jo Jo Jo 

= [ T dT M P (T - TV) j io T dxj ^ 9 ( X ~ T ^ 
Jo 

Jo 

&T X f(T 2 - Ti) e'fo dxj(x)g{ X - Tl ) f ^ ^ ^ 

Jo 

where we have used the identity jfi <*> = ftf =1 e * / T <^'M <^-t ; ). Note thftt . f 
j(x) = 0, the functional Zj^\j = 0] reduces to Pa/(T), as stated above. This equation can 
be written in a more condensed form if one introduces the following set of functionals, 
defined recursively by 

lo(j,u) =/(«) 

WU, «) = f dv f(u - v) e l Io «M»)*<«-«) V ) . (28) 

Jo 

Using this notation, it is then easy to see that one can write the functionals Zj^[j] in 
the form 

Z \j] =P {T) 

ZjM = / T d u Po{T - u) e'fi **m 9 {*-u) z ^ u) if AT > . (29) 
Jo 
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Since we wish to compute the functional Z\j\, we need to performed the summation 
defined in equation (|24|) . If one now introduces the functional 



l(j,u)= E WU,u), (30) 

Af=0 

then one can show from equation (|2*Hj) that l(j,u) obeys the integral equation 

l(j )U ) = f( u ) + f dvf(u - v) e*So tej^gix-v) ^- ^ _ ( 31 ) 



Now, using the recursion relation ()29|). it is easy to show that Z[j] can be expressed in 
terms of u) by 

ZU] = P Q (T) + f duP (T - U) e^o dxj{x)g{x-u) l( j ■ s _ (32) 

Jo 

As they stand, equations (|3*T|) and (|3^|) are of little use, since it is usually not possible 
to solve (f3~TJ) for l(j,u). However, if one takes j(x) = in equation (J3T|) . one obtains, 
after changing the integration variable from i> — > u — v, the equation 

U0 u ) = f(u) + [ U dvl(0, u-v) f(v) , (33) 
Jo 

with the initial condition Z(0, 0) = /(0). But this is precisely equation (JTHJ) for the 
renewal function l(u), with the same initial condition. Therefore, we conclude that 
1(0, u) = l(u), a result which justifies the notation we are using. Now, we express f(u) 
and f(u — v) in equation (|3*T|) in terms oil(u), using (}3*3*j) . After interchanging the limits 
of integration in the resulting double integral and using (J2J, one obtains the integral 
equation 

l(j, U ) = 1(U) + £ dvl(u - V) (jto dxj(x)g(x-v ) _ ^ ^ ? ( g 4) 

which shows explicitly that l(j,u) reduces to l(u) when j(x) = 0. Now, substituting 
Po(T) = 1—Jq du f(u), Pq(T—u) = 1— /J dv f(v—u) inequation and interchanging 
the limits of integration in the resulting double integral, one obtains, after using equation 

m 

Z[j] = 1 + £ du (V/o T »(*-«) _ ; (35) 

which shows explicitly that Z[j — 0] = 1. Equations (}3"4"j) and (|33j) are the main results 
of this section. In fact, these two equations are equivalent to the expansion of the 
probability-generating functional Z\j] in terms of the conditional probability densities 
for a renewal process, see [HO] for details. The above form is more convenient for explicit 
calculations, as shown below. 

Before we discuss the general solution of (J34|) . let us apply our results to compute 
the functional Z[j] for the particular case of the simple Poisson process. Given that 
l(u) = 1/r, one concludes, substituting this result in equation (jMj), that l(j,u) obeys 
the integral equation 

l(j, U) = -j- + £ y ( e'fo dxj(x)g( X -v) _ A v) (36) 
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Differentiating this equation with respect to u, one obtains the first-order linear 
differential equation 

^jA = I (jfi *>i(»)*(»-«0 _ i M . (37) 

with the initial condition 0) = 1/r. Such an equation can be quickly integrated by 
separation of variables, yielding the solution 

l(j,u) = lexp ( [ U *L ( e <£***W"» -l)). (38) 



10 T 

This solution, when substituted in equation (J3*5|) . yields 

Z[j] = exp ^ ^ (V Jo" - l) ) , (39) 

which is a well known result, being given in, e.g. [3*T] . However, this is in fact the 
only case for which a closed expression for Z[j] can be found, as one can no longer 
reduce the integral equation (|3*4*)l to a first-order differential equation if one chooses a 
different waiting-time distribution. Neither can one solve ()34|) using a Laplace transform 
because the source term e jo j{x) g{x-v) j n j g no ^ time-translation invariant. Such 
difficulties stem from the fact that equation (J3*4"|) is a Volterra integral equation for which 
no analytical solution is known in the general case. 

Nevertheless, one can show that all correlation functions can be obtained from 
equations (J3*4*)l and (J33j) in a closed form, provided one knows the renewal function l(u). 
This point can be made clearer through the discussion of two examples which we will 
use extensively in our applications of renewal processes to the study of the dynamics 
of molecular motors. We wish to compute the average value {J{ti)) and the two-point 
correlation function {J{t\) J^))- In the first case, one has 



3=0 




du g{h - u) jfi l(j, u ) 

( J Jo dxj{x)g(x~u) _ A 

V J 61(h) 



3=0 



T 

dug{t x - u)l(u) , (40) 

o 

since the second term is zero if j(x) = and l(j = 0,u) = l(u). Therefore, the average 
value (J(ti)) can be expressed as an integral involving the known functions g(u) and 
l(u). In particular, if J(t\) = N(ti), which is the number of increments which has 
occurred up to time t± or the displacement of the particle up to time t\, then one has 
g(u) = 6{u) and one obtains 

(N(ti))= [ H dul{u) = m(ti), (41) 
Jo 

which is just the definition of the renewal function m{t\) given in the previous section. 
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The procedure used to compute the two point correlator (Jfa) J fa)) is analogous, 

^tional 

s 2 m 



but it also involves the functional derivative f l M' u \ ■ One has 



{J fa) J fa)) 



T 

dugfa — u) gfa — u) l{u 

T 



Sjfa)5j(t 

i / dugfa — u) j —if dug{t 

JO ditto) , „ J0 



6l(j,u) 

2 - U 



3=0 



,(42) 

3=0 



since the term which involves the functional derivative sS^jr^) * s zero ^ •^* E ) = ^' 
Now, in order to compute the functional derivative ^"j* at = 0, one functionally 
differentiates equation (J54*|) with respect to jfa) and sets = 0. One obtains 

l {[ dvl ( u ~ v ) 9 fa ~ v) e^o **j(x)g(*-v) 

+ r dvi(u - v) (e*i: - 1) ^4 1 

i /* cfa Z(« - v) gfa - v) l(v) , (43) 



Sjfa) 



j=0 



since the second term is again zero if jfa) = and l(j — 0, v) = l{v). The expression 
for the derivative j$0; is simply obtained from equation ()43j) by interchanging t% by £2- 

It should be clear from equation (}4T)j) that in order to compute the m-th variational 
derivative of Z[j] at jfa = 0, one needs only the derivative of order m — 1 of l(j,u), 
because the term which involves the m-th derivative of l(j, u) is zero in that case, since 
e ij dxj(x) g (x-u) _ i = if = 0. More importantly, it should also follow from the 
structure of equation (|43|) that in order to compute the m-th variational derivative of 
l(j,u) at jfa) = 0, one only needs the derivative of order m — 1 of l(j,u), because the 
term which involves the m-th derivative of l(j,u) is zero for the same reason. This 
implies that one can express such variational derivatives in terms of quantities which 
were previously calculated, assuming that we know the function Ifa). Thus, the solution 
of equation (|T9| is sufficient to compute all correlation functions in a recursive manner, 
sec [30| for a rigorous proof. 

Finally, one can obtain a closed expression for {Jfa) Jfa)) by substituting the 
result for f^'^ given in equation (gSJ) in equation (pEZj) . One has that (Jfa) Jfa)) 
is given by 



(Jfa) Jfa)) = ^ du gfa - u) gfa - u) l(u) 



+ 



/ dugfa — u) \ dvl(u — v) gfa — v)l(v) 
Jo Jo 



T 



+ dugfa- u) l dvl(u - v) gfa -v)l(v) . (44) 
Jo Jo 

If Jfa) = Nfa), Jfa) = Nfa), then gfa) = 6(u) and we obtain for the correlation 
function (Nfa) Nfa)), the result 

/•min(ti,ta) pti i-minju^) 

(Nfa) Nfa)) = / dul(u) + du dvl(u-v)l(v) 
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+ du dv l(u - v) l(v) . (45) 



In the context of the study of molecular motors, it is more appropriate to consider 
the mean-square deviation of the displacement of the particle between time t 2 and ti 
(ti > t 2 ), i.e. the correlation function ((iV(ti) - iV(t 2 )) 2 ) C onn = ((JV(*i) - N(t 2 )) 2 } - 
(N(ti) — N(t 2 )) 2 , as explained in sectional If we use equations (j4~Tj) and (J43j) . we obtain, 
after some trivial manipulations, 

((N(ti) - N(t 2 )) 2 ) conn = f 1 dul{u) + 2 t du l(u) r dv {l{v -u)- l{v)). (46) 

Jt2 Jt2 Ju 

In the case of a Poisson process, l{u) = 1/r and this expression reduces to the first 
term, i.e. we obtain ((N(ti) — iV(t2)) 2 )conn = (£i — ti)l T i which is a well-known result 
for the simple Poisson process. The expression for the full functional obtained above in 
equation (JH9*j) allows one to compute correlation functions of any order, but only in this 
simple case. 

A particular important limit of (|46|) is the one in which t 2 ,ti — > oo, but in which 
the time difference t — t\ — 1 2 is kept finite, since C(t) = lim tl t2 ^ 00 ((A^(t 1 ) — iV(t 2 )) 2 ) C onn 
becomes a function of t only, i.e. one recovers a form of time translation invariance. 
In this limit, l(u) — > -tV and one obtains, substituting this result above, the following 

result for C(t) 

CW = ^ + ^jf*.(*-«)(l(«)-^). (47) 

If t <C (r), it is easy to see that the second term of this equation is at least 0(t 2 ) and 
therefore C(t) m -jK-, if t <C (r). This observation justifies the first line of equation (J3J). 

In order to extract the long-time behaviour of C(t), one needs to consider instead the 
Laplace transform C(s) of C(t). One has, from equation (jJJJ), 

C( S) 1 (lip *), (48, 

where we have used equation (J2~U|) to express /(s) in terms of f(s). If /(s) is analytic 
a s = 0, i.e. if all moments of the distribution /(£) exist, one can write, for small s 

/ r 2\ / r 3\ 

/(s) = 1 — (t) s + -^p s 2 — ^-gp- s 3 + . . .. Performing a Taylor expansion of at s = 
using such a result for /(s) and keeping only the most divergent terms, one obtains 

C(s) « -^— 2 + - , (49) 

with r = ij2 )~^ and C = - ^X, where (r 2 ) = / °° d(( 2 f(() and (r 3 ) = 

(r) 2(r) s(r) 

Jo°° ^CC 3 /(0- Furthermore, one can write C in the form given by equation (jSJ), using 
the definitions of the connected second- and third-moment r and S, given above. The 
behaviour of C(t) for t ^> (r) is determined by the behaviour of the Laplace transform 
C(s) at small s, as given by (fUJj) and one can directly identify the coefficients of the 
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most divergent terms of C(t) from those in equation (jjHJ) . One thus concludes that C(t) 
is given, in this limit, by the second line of equation (0J). 

Let us now consider the case in which the waiting-time distribution of the renewal 
process is given by a convolution of M. Poisson processes, each of which is supposed to 
represent a rate-limiting step of the molecular motor. Mathematically, this is expressed 
by saying that such function, which we denote by /.m(C), is given recursively by 

f M (0 = f — f M -i(v) ^ M > 1 , (50) 

JO T M 

where f'M-iiv) is the waiting-time distribution of a process with Ai — 1 rate-limiting 
steps and /i(C) — e^^ Tl jr\. Given that the Laplace transform of a single Poisson process 
is given by fi{s) = (1 + Tis)^ 1 , one concludes from equation (|5T)|) and the convolution 
theorem that the Laplace transform of fjuiO is given by 

M i 

f~M(s)=l[-—. (51) 

Since (r) = —f'(0), (r 2 ) = f"(0) and (r 3 ) = — /^(0), one concludes, differentiating 
equation (joTj) . that (r), r and C are given, respectively, by equations (JTJ), (jSJ) and (0) 

with 5 = . j^ 1 1 a . The form which /^(s) takes for = 2, 3 even allows one to 

(E,=i r 

invert the Laplace transform l(s) and obtain l(u) in these two cases. One can therefore 
determine C(t) for arbitrary time t, using equation (|47|). This calculation is explicitly 
performed in |Appendix as stated above. 

In the above derivation, we have assumed that all moments of the distribution f(Q 
exist, in other words that the Laplace transform f(s) is an analytic function at s = 0. 
Such assumption is enterily justified when discussing the motion of a molecular motor 
in a homogeneous environment, as is the case of myosin V or kinesin. On the other 
hand, if one considers the motion of tracer particles in a rapidly rotating fluid jHI], one 
is led to consider a renewal process with a waiting-time distribution with fat tails, which 
reflects the diverging sticking-times of the tracer particles. In a biological context, some 
of these results may also be applicable to the study of the motion of RNA-polymerase 
along a DNA-strand [Zj. 

It is beyond our means to provide a complete discussion of this issue, However, one 
may consider a very simple model for a waiting-time distribution with fat tails, namely 
we take f(() = 7^^ v , where v > 1, such that /(£) is normalisable. For 2 < v < 3, 
the assymptotic decay of f(() is the same as that of a stable Levy distribution The 
Laplace transform of /(£) can be readily calculated, at least in the limit of small s, and 
one may carry through a significant number of calculations with this simple model. 

One should start by considering the consequences of taking v < 2, in which case 
the distribution f(() has an infinite first moment (r). In such a case, one can neither 
obtain equation flj7j) from equation nor equation (JHJ) from equation (see below), 
since we are assuming that (r) is finite in the derivation of these results. Likewise, the 
equations that follow from equation (}T7j) . such as equation ([48)1. and from equation (0), 
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such as equation (J7J), are equally invalid. In physical terms, the assumption that fails is 
that one can consider the existence of a time-translation invariant regime, ti,t 2 ^> (r), 
for which the correlation functions ((N(ti) — N (t 2 )) 2 ) conn and S(q,ti,t 2 ) depend only 
on t = t\ — t 2 , since (r) is formally infinite. 

In the case in which v > 2, a time-translational regime does arise and our results 
hold through, except equation (|4"9j) which relies on the assumption that f(s) is analytical 
at s = 0. In particular, one may substitute the Laplace transform f(s) in equation (|47|). 
and extract the asymptotic behaviour of C(t) at large t, by considering the limit of small 
s. The results we have obtained are simply a limiting case of those obtained by Weeks 
et al. [SI], and we will therefore only quote the final result. 

One can show that at large t, C(t) has the following asymptotic behaviour 

( t 4 ~ u if 2 < v < 3 
C(t)~| tint iff = 3 , (52) 

t^t t if v > 3 
I (r) 

with r = > 1 being the motor's (finite) randomness coefficient if v > 3. 

The result obtained in the parameter region 2 < v < 3, in which the second 
moment of the distribution, or r, is formally infinite, can be interpreted as a super- 
diffusive behaviour of the tracer particle. Such a behaviour is also present in the model 
studied by Kafri and co-workers [7| 

We now consider the derivation of the Laplace transform of the density-density 
correlation function S(q,ti,t 2 ) = ( e~ %qd<yNtytl ^ Nt - t2 ^ ) in the long time limit £1,^2 — > 00, 
which is given in equation (JHJ). Such a correlation function can be obtained from the 
general expression for the probability-generating functional Z[j] by a judicious choice 
of g{u) and of the source function j(x) in equations ()22I24|) . If we take g{u) = 6{u) and 

j(x)=qd(5(x-t 2 )-5(x-t 1 )) (53) 

then it is easy to see that Z[j] reduces to S(q,ti,t 2 ) for this particular choice of the 
functions g{u) and j(x). It also follows from equation (|3~K|) that S(q,ti,t 2 ) obeys the 
following integral equation 

S(q, t x , t 2 ) = l + (e~^ - 1) f d( l(q, t 2 + (), (54) 

Jo 

where t = t\—t 2 and where l(q, t 2 + () is given by the solution of the integral equation 

Kq,h+C) = l(t 2 +C) + (e- iqd -l) f drjl(C-v) Kq,t 2 +v) if < C < *-(55) 

J 

This equation follows from ()34)1 for the particular choice we made for g(u) and j(x). Such 
a system of equations can be solved using Laplace transformation. Following Godreche 
and Luck (2Z|, we define the following Laplace transforms 

l(s;t 2 ) = / d(e- s U(t 2 + (), (56) 



l(q,s;t 2 ) = / dCe-<l(q,t 2 + C), (57) 



00 



S(q,s;t 2 )= / dte~ st S(q,t 2 + t,t 2 ) ■ (58) 
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Substituting equations (|54I55|) in (|57I58|) . one obtains, after some trivial manipulations, 
the following result for S(q, s; t 2 ) 

1 / e~ iqd - 1 ~ \ 

% , S;t2 ) = -(l + i + (i _ e _ w)f(s) /( S;fa )j (59) 

where l(s) is the Laplace transform of l(u), given by ([20)1 . In the long time limit t 2 — > oo, 
Z(i 2 + C) "~ > T^y an d l(s]t 2 ) = which is independent of £2- Therefore, S(q,s;t 2 ) 

is also independent of t 2 and we write it simply as S(q, s) in this limit. Substituting 
l(s) by its expression in terms of f(s), as given by equation (J2*0"|) . in equation (J59|) . we 
conclude that S(q, s) is given by equation (jUJ) 

We will now consider the relation between the Fourier and the Laplace transforms 
of S(q,t). From its definition, one immediately concludes that for t > 0, 



S(q,-t)=S(-q,t) = S(q,t). (60) 
If we substitute this identity in the definition of S(q,u), we obtain 

dte iujt S(q,t) + / dte twt S(q,t) 
-00 J 

= / dte^ S(q,t) + / dte iujt S(q,t) 
Jo Jo 

= 2ReS{q,s = -iu) , (61) 

where we have performed the change of variable t — ■> — t in the first term of the first line 
of (}6T|) . Equation ((7j) then follows from (jUJ) and from (JBTj) . Finally, one obtains (jSJ) by 
substituting f(—iu) = (1 — zcjt) -1 , valid for a simple Poisson process, in ((7|). 

Before closing this section, we will briefly explain how one can obtain the expression 
for the spin-spin correlation function (a tl ■ <?t 2 ) in the time domain in the cases in which 
the waiting-time distribution of the renewal process is given either by a single Poisson 
process, as in (|TT|) . or in which the waiting time-distribution is given by the convolution 
of two Poisson processes with the same characteristic time r/2, as in (|12jl . leaving 
the details to Appendix C If the waiting-time distribution f(() is a simple exponential, 
then, substituting the result for its Laplace transform f(s) in equation (flT)|) . with Q = 3, 
one observes that the Laplace transform of the spin-spin correlation function has a simple 
pole, and one can directly read (a tl ■ a t , 2 ) from it. If the waiting-time distribution is 
given by the convolution of two Poisson processes with the same characteristic time, the 
Laplace transform of (a tl ■ a t2 ), as given by equation (JTUJ), with Q = 3, has two simple 
poles and this expression can always be written in terms of partial fractions involving 
one or the other of these poles. In that case, one can also read (a tl ■ <?t 2 ) directly from 
it. 



5. Conclusion and outlook 



We have modelled the dynamics of a processive or rotary molecular motor as a renewal 
process, in line with the work of Svoboda and co-wokers. Using a functional technique, 
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we have computed the mean-square deviation of the distance travelled by a processive 
motor and extracted its asymptotic limit at large times. For renewal processes composed 
of two or three Poisson substeps, we have computed this function outside the asymptotic 
regime, given the relevance of such processes for the study of myosin V and kinesin. It 
follows from our results that the measurement of this correlation function would permit 
one to extract additional information concerning the time constants which characterise 
the motor's chemical cycle. 

We have also used the same functional method to compute the density-density 
correlation function of an ensemble of independent processive motors, which can be 
measured using the experimental techniques developed by Cappello and co-workers. 
We have also shown that in a particular limit such a function reduces to the spin-spin 
correlation function of a 'random-clock' model that has applications to the dynamics 
of rotary motors, such as Fl-ATPase. The measurement of this correlation function, 
followed by its Fourier transformation to frequency space, would permit one, in both 
cases, to fit such a quantity to the theoretical results discussed above. Thus, one could 
obtain all the time constants characterising the molecular motor's chemical cycle, even 
in the case of chemical cycles composed of a large number of rate-limiting substeps. 

As possible avenues of future research, one can indicate at least two experimental 
issues that still need to be addressed. The first experimental issue is the repetition of 
the experiments of Cappello and co-workers with an interference mask with a smaller 
period or the use of myosin V, rather than kinesin, as the subject of study of such 
an experiment, which would allow the measurement of the density- density correlation 
function away from the limit of long wavelengths, where the motor simply behaves 
as a Brownian particle. The second issue is the direct measurement of the spin-spin 
correlation function in Fl-ATPase using fluorescence microscopy and how one can 
improve the time resolution of the present technique. 

In closing, we may say that the results presented in this paper can be used, in the 
different experimental contexts to which they apply, to provide for a precise fitting of 
the time constants associated with the rate-limiting steps of a molecular motor chemical 
cycle. It remains to be seen to what extent such information is of crucial importance 
to the understanding of the chemical kinetics of processive or rotary molecular motors, 
or if a qualitative understanding of the nature of the chemical cycle is by itself sufficient. 
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Glossary 



Bead-motor assay: Experimental apparatus where a spherical plastic bead of 
micrometric size is coupled to a single molecular motor. The position of the bead 
can be controlled with high precision using two intense laser beams focused at a spot 
(known as optical tweezers). Due to its optical properties, the bead is attracted to the 
focus of the two beams and may be controlled at a distance. 

Molecular motor stepping: Directed motion of a processive molecular motor (e.g. 
myosin V, kinesin, dynein) along a specific molecular track in the cell (e.g. actin in 
the case of myosin V, tubulin in the case of kinesin and dynein). Such a motion is 
composed of individual steps of fixed length, which occur at random times, determined 
by the chemistry of the process. 

Probability-generating functional: Mathematical object which encodes in itself 
all the information that can be obtained (measured) from a random process. In the 
main text, we have computed the explicit form of such an object for a renewal 
process where the interval distribution function was a simple exponential. From such 
an expression, one can obtain all the correlation function pertaining to such a renewal 
process. Thus, the ultimate goal of applied probability theory is the computation of the 
probability-generating functional for a given random process. All too often, one has to 
content oneself with well less than that. 

Renewal process: Random counting process where an integer variable is increased 
by one unit at random times, the statistical distribution of the length of time intervals 
between sucessive counting events being a known function. 

Appendix A. Mean-square deviation of the displacement for a renewal 
process composed of two or three Poisson substeps 

In this appendix, we derive the explicit form, valid for arbitrary time t, of the mean- 
square deviation of the distance travelled by the motor C(t), for a renewal process whose 
waiting-time distribution is given by the convolution of two or three Poisson processes. 
As can been seen from equation (|47|) . the knowledge of the renewal function l(u) suffices 
to determine C(t). 

We start with the simplest case, namely the case in which the waiting-time 
distribution is given by the convolution of two Poisson processes. In that case, the 
Laplace transform of the waiting-time distribution is given by f(s) = ^ 1+ST1 ^ l+ST2 ^ , as 
follows from equation (foTj). Substituting this result in equation (j2Uj). one can show that 
the Laplace transform of the renewal function can be written in terms of partial fractions 
as 



From this equation, one can directly determine the inverse transform l(u), which is given 
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by 



l(u) 



-I ± + ± )u 



n + t 2 



(A.2) 



which tends to l{u) 



in the limit of large u, in agreement with Blackwell's 



renewal theorem. Substituting this result for l(u) in equation (j4*Tj) and performing the 
integration over u, one obtains for C(t) the result 

rt ItItI 



C(t) 



where r = ^+^2 is the randomness parameter for a motor whose chemical cycle is 
composed of two rate-limiting steps. It can be easily checked that this function has the 
correct asymptotic forms at small and large t, has given by We have plotted this 
function in figure ^ with n = T2 = 1/2, i.e. when r = 1/2. 

In the case of a renewal processes whose waiting-time distribution is given 
by the convolution of three Poisson processes, it follows from (JoT|) that f(s) = 
rn — rrrr — rrn — r- Substituting this result in equation (I2U1) . one can show that the 
Laplace transform of the renewal function can be written as 



+ 



T1 + T2 (ti + t 2 



1 



-(-+-)t 



(A.3) 



1(8) 

where 7± = 



1 



n+r 2 + r 3 



7+ 



7- 



7+ -7- 



s + 7- 



7+ - 7- 



4 



T1T2 Tlf"3 T 2T3 



(A.4) 
One can also 



f^ + i + ^)± 

V n t 2 r-i J 

directly determine the inverse transform l(u) from this equation, which is given by 

1 ( 7+ , A „- y _ u) 7- 



l(u) = - 

Tl+T 2 + T 3 

This function tends to l{u) 



-7+ u\ 



,7+ 



. (A.5) 

7- 7+ - 7- J 

— — —. — in the limit of large u, which also agrees with 
T1+T2+T3 

1 7 ix \v< 117 renewal theorem. Substituting this result for l{u) in equation (|4*T|) and 
performing the integration over u, one obtains for C(t) the result 




- 7 _t 



where r 



(A.6) 



is the randomness parameter for a motor whose chemical cycle is 



(T1+T2+T3) 2 

composed of three rate-limiting steps. It can be easily checked that this function has 
the correct asymptotic forms at small and large t, has given by (jlj. We have plotted 
this function in figure |2l with T\ = r 2 = r 3 = 1/3, i.e. when r = 1/3. 



Appendix B. Density-density correlation function in Fourier space for a 
renewal process composed of two or three Poisson substeps 



All one has to do is to substitute the appropriate form for f(—iuj) in equation (J7J). For 
a renewal process whose waiting-time distribution is given by the convolution of two 
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Poisson processes, f{-iu) = {1 _ iuJTl ) (1 _ iuJT2) and S{q,u) is given by 

2(n + r 2 ) (1 - cos(qd)) (r + o^ 2 ) ^ 
ig ' J (cu(r 1 + r 2 )-sinM)) 2 + (l-cosM)-r 1 r 2 ^) 2 ' l ' y 

r 2 +r 2 - 2 - 2 

where r = 7^- — ^ and a 



(T1+T2) 2 (1-1+T2) 2 ■ 

For a renewal process whose waiting-time distribution is given by the convolution 

of three Poisson processes, f(—iu) = ( 1 _^ Tl )( 1 _ t 1 a;r2 )( 1 _^ T - 3 ) an d S(q,u) is given by 

S(q, u) = 2(n + r 2 + r 3 ) (1 - cos(gd)) (r + oxu 2 + /fo 4 )/iV (B.2) 

where 

iV = [u(ti + r 2 + r 3 ) - sin(gd) - rir 2 r 3 u; 3 ] 2 

+ [1 - cos(grf) - (nr 2 + rir 3 + r 2 r 3 ) w 2 ] 2 , (B.3) 

^-2 1 ^-2 j ^-2 ^_2^_2_|_ 7 _2^_2_|_^_2^_2 ^_2^_2^_2 

aIld r = (n+ra'+ra 3 ) 2 ' a = 1 (r 1+ r 2 |r 3 ) 2 " aIld ^ = (r 1 +r 2 Vr 3 ) 2 " 



Appendix C. Spin-spin correlation function for a renewal process composed 
of two Poisson substeps 

We first obtain the spin-spin correlation function in the time domain, as given by 
equation (fT2|l. if the two time constants are equal, i.e. if T\ = r 2 = r/2. In such a 
case, the Laplace transform of the waiting-time distribution is given by f(s) = / 2 ) 2 • 
Substituting such a result in equation (fTUj) . one can write the Laplace transform of the 
spin-spin correlation function in terms of partial fractions as 



C {a h ■ ot a ) = Re 



3 1 11 



+ - 



(G.l) 



4 s + \ (1 - e"* 71 "/ 3 ) 4s + ^(l + e-"/ 3 ) 

From such a result, one can immediatelly read its inverse Laplace transform, which is 
given by equation (JT2J) ■ Note that if the renewal process is a simple Poisson process, 
f( s ) = 14^7 an d the Laplace transform of the spin-spin correlation function involves a 
simple pole, which make its inversion trivial and we obtain ([lip. 

For a general choice of the time constants T\ and r 2 , one can also write the 
Laplace transform of the spin-spin correlation function in terms of partial fractions, 
which permits to invert such a transform. However, extracting the real part of such an 
expression becomes a cumbersome exercise, albeit a trivial one. It is simpler, also from 
the viewpoint of the experimental fitting of the data, to consider instead the Fourier 
transform of the spin-spin correlation function, which is defined as FT(at 1 ■ &t 2 ) = 
J"_ro dte iuJt (a tl ■ <?i 2 ), with t = ti — t 2 . Since one can write the spin-spin correlation 
function as (a tl -Ot 2 ) = |( S(2ir/3, t) +S(2ir/3, —t) ) (with d = 1), the Fourier transform 
of this function is given by 

FT(a h -P t2 ) = \ (S(2tt/3, u) + S(2tt/3, -u)) . (C.2) 
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where S(q = 2tt/3,uj) is given by ([B.ip . valid for a renewal process composed of two 
Poisson substeps. We finally obtain 



FT(a tl -a t2 ) = -(r l + r 2 )(r + auj 2 ) 



1 



+ 



(w(n + r 2 ) - v/3/2) 2 + (3/2 - nr 2 uj 2 ) 2 
1 



(C.3) 



where r 



and a 



w(ri + r 2 ) + v/3/2) 2 + (3/2 - nr 2 cu 2 ) 2 

_2_2 

. We plotted this function for three different choices 



(T1+T2) 2 " Jli "- (ri+T2) 2 

of the time constants Ti and r 2 , such that T\ + r 2 = 1, in figure |U 
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